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1. Introduction 

Coupled nonlinear system may exhibit a remarkable phenomenon called chaos 
synchronization [H |2]. The individual systems synchronize to a common chaotic 
trajectory. This phenomenon has received much attention, due to its potential 
applications in secure communication P, H]. For technological applications semi- 
conductor lasers are promising systems to implement secure communication schemes 
using chaos synchronization, because they exhibit fast dynamics, they are cheap and 
one could utilize the existing telecommunication infrastructure for these lasers |5]. 

However, due to the fast dynamics of the lasers, propagation distances of already 
a few meters introduce non-negligible delay times in the coupling. The synchronization 
of delay-coupled systems in general [HI El El E] and in particular delay-coupled lasers 
[TOl im [T2I [T3l [HI [151 UHl [13, [IE] has thus been a focus of research in nonlinear dynamics 
during the last decades. 

When the lasers are coupled ail-optically, not only delay effects are important, 
but also the optical coupling phases of the coherently coupled electric fields play an 
important role. Coherent coupling may result in constructive or destructive interference 
of incoming signals. When the lasers are synchronized, this interference can occur even 
if the coupling distance is much larger than the coherence length of the beams |19) . 

It has been demonstrated experimentally and numerically that by tuning coupling 
phases one can adjust the level of synchronization ranging from perfect synchronization 
to almost no correlation [20] • However, the impact of these coupling phases and the 
resulting interference conditions have so far not been thoroughly investigated. 

In this work we derive and discuss synchronization conditions for ail-optically 
coupled lasers. These conditions ensure the existence and stability of a completely 
synchronized solution, i.e., perfect synchronization of identical lasers. We confirm the 
analytical results by numerical simulations. For our stability analysis we employ recent 
results concerning the master stability function (MSF) [21j in the limit of large delay 
times 1221 [23] • 

The paper is organized as follows. In Sec. [2]we briefly review previous results about 
synchronization in delay-coupled networks in a general context. These general results are 
then applied to ail-optically coupled lasers. In particular, we discuss the case of two ail- 
optically coupled lasers in detail in Sec. [3j Here, we consider different coupling schemes 
(unidirectional and bidirectional, open- and closed-loop) and derive the corresponding 
synchronization conditions. We then discuss synchronization in larger laser networks in 
Sec. |4] focusing on an experimentally feasible setup. Finally, we conclude with Sec. [5j 

2. Synchronization of delay-coupled systems 

Consider a delay-coupled network of identical units 



N 




(1) 
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where Xk G is the state vector of the k-th node {k = 1,...,N). Here, / is a 
function describing the local dynamics of an isolated node, Gkj is a coupling matrix 
that determines the coupling topology and the strength of each link in the network, h is 
a coupling function, and r is the delay time in the connection, which is assumed to be 
equal for all links. A necessary condition for perfect synchronization is that the matrix 
G has a constant row sum 

N 

CI: a = J2Gk,, (2) 

i=i 

independent of k. This condition ensures that an invariant synchronization manifold 
exists. Note that in this case the row sum cr is an eigenvalue of the matrix G 
corresponding to the eigenvector (1, !,...,!) of synchronized dynamics. For a given 
constant row sum a the dynamics in the synchronization manifold of Eq. ([T]) is 

^(t) = f[x{t)] + ah[x{t - r)] . (3) 

where x{t) = Xi{t) = Xi{t) = . . . = x^it) denotes the synchronized solution. 

Note that condition CI addresses mathematically perfect synchronization. For 
realistic systems it is important to investigate what happens to the synchronization 
manifold under small perturbations of the perfect setup, such as parameter mismatch 
of the individual elements or coupling strength mismatch leading to a slightly broken 
condition CI. On one hand, there has been some recent progress in addressing parameter 
mismatch within the MSF framework fl^ [2S] and it might be interesting to apply these 
new methods to the laser system. On the other hand, slight perturbations of the perfect 
setup are also related to the question whether the synchronization manifold is normally 
hyperbolic Although these questions are important for realistic systems, we will 

focus our analytic investigations on the case of perfect synchronization to highlight our 
main message: the importance of the coupling phases. 

As we will discuss below, the necessary condition CI becomes quite involved for 
the case of optically coupled lasers, because here the matrix G is complex valued due 
to optical phase factors in the coupling and, additionally, the lasers can have relative 
phase shifts, which effectively allow the system to adjust the coupling matrix to some 
extent. 

If condition CI is satisfied and a synchronization manifold exists, the stability 
problem of the synchronized solution can be approached using the MSF |21j . The MSF 
depends on a complex parameter re^"^ and is defined as the largest Lyapunov exponent 
A(re*^) arising from the variational equation 

m = Df[x{t)] m + re'^Dh[x{t - r)] i{t - r) , (4) 

where x{t) is the synchronized trajectory of the system determined by Eq. ^ and Df 
and Dh are Jacobians. 

The synchronized state is stable for a given coupling topology if the MSF is 
negative at all transversal eigenvalues 7^ of the coupling matrix (A(7fc) < 0). Here, 
transversal eigenvalue refers to all eigenvalues except for the eigenvalue a associated 
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to perturbations within the synchronization manifold with corresponding eigenvector 
(1, 1, 1). 

The dynamical time scale of a semiconductor laser is given by the relaxation 
oscillation period Tro- Typical values of Tro for semiconductor lasers are between 
lO'^^s and 10~^s. This time scale corresponds to about 3cm - 30cm of optical path- 
length in air. The coupling distances between the lasers is usually of the order of meters 
or even much larger [5]. In this case the coupling delay time is much larger than the 
intrinsic time scale r ^ Tro of the lasers. The case of large delay times is an import 
hmit for delayed systems [271 |28l [291 |30l [3ll [32] . 

We recently showed [25] that in any network the stability problem for the 
synchronized state is drastically simplified in this limit of large delay times: 

• The MSF is rotationally symmetric around the origin in the complex plane, i.e., 
A(re*'^) is independent of ip. 

• If A(0) > 0, then X{re'^) = A(0) for all r (and ^j). 

• If A(0) < 0, then the MSF is monotonically growing with respect to the parameter 
r and there is a critical radius tq, where it changes sign (A(ro) = 0). 

Recently, this structure of the MSF has been been confirmed experimentally [33] and 
even utilized to predict synchronizability of a network from a simpler motif [331 

Note that the two cases A(0) > and A(0) < have recently been discussed in 
more detail [23]. The case A(0) < is called weak chaos, since in this case the maximum 
Lyapunov exponent scales as ~ 1/r for r — > oo. The case A(0) > 0, on the other hand, 
is called strong chaos since the maximum Lyapunov exponent scales as 0(1) for r — )■ oo. 
It has been shown that for large delay networks can only exhibit chaos synchronization 
in the regime of weak chaos [22^ [23] . 

The structure of the MSF allows us to draw general conclusions about the 
synchronizability of a given network topology. In particular, chaos synchronization can 
only be stable if a is the eigenvalue of G with largest magnitude, i.e., the transversal 
eigenvalues 71, 72, . . . 7Ar_i have smaller magnitude 

C2 : |7„| < \a\ . (5) 

This condition (C2) is necessary for chaos synchronization (|cr| > Tq) and it is sufficient 
for synchronization on a periodic orbit {\a\ < tq) [22]. In fact the smaller the magnitude 
of the transversal eigenvalues, the easier it is to synchronize the system. 

For given synchronized dynamics (Eq. ([3])), i.e., given system parameters and row 
sum a, one can calculate numerically the critical radius Tq, which then provides a 
necessary and sufficient condition (C3) for synchronization (provided CI is fulfilled) 

C3 : |7„| < ro . (6) 

While condition C3 is necessary and sufficient, its disadvantage is that one needs to 
calculate tq explicitly for the particular synchronized dynamics. In contrast, condition 
C2 is only necessary for chaotic synchronization, however, its strengths lies in the fact 
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Figure 1. Schematic coupling scheme of two delay-coupled lasers with delayed self- 
feedback. Each connection has a coupling strengths Kkj and a coupling phase (jj^j- We 
consider the case when all delay times are equal. 



that it depends solely on the coupling topology, i.e, the eigenvalues of G, and not on 
the particular dynamics of the system. 

The aim of this paper is to apply the synchronization conditions C1-C3 to a system 
of two delay-coupled lasers and to laser networks and discuss the consequences for chaos 
synchronization. In particular, we show how condition CI gives complicated conditions 
for the coupling phases of the lasers. 



3. Two Lasers 



In this section we consider two semiconductor lasers that are delay-coupled to each other 
with a coupling delay and additionally receive self-feedback with the same delay time 
r. The basic coupling scheme is depicted in Fig. [T| The coupled system is described by 
dimensionless rate equations of Lang-Kobayashi type 

El = 1(1 + ia) [G{ni, El) - 1] El 

+ Kiie'^'^Ei{t - r) + Ki2e''^'^E2{t - r) , 
E2 =\{l + ia)[G{n2, E2)-l]E2 

+ K22e'^''E2{t - r) + K2ie''^''Ei{t - r) , 
Till =p-ni- G{ni, Ei)\Ei\^ , 

Th2=p-n2-G{n2,E2)\E2\\ (7) 

where Ek and Uk are the normalized complex electric field amplitude and the rescaled 
inversion of the k-th laser, respectively, and a is the linewidth enhancement factor, p 
is the normalized pump current in excess of the threshold, and Hij, (pij are the coupling 
amplitudes and phases, respectively, as shown in Fig.. [TJ The gain is modeled by 

that takes into account gain saturation effects. Throughout this paper we choose the 
following model parameters for our numerical simulations, unless stated otherwise: Ratio 
between carrier and photon lifetime T = 1000, p = 0.1, a = 4, gain saturation /i = 0.26. 

The optical coupling phases (pij are determined by the optical path lengths of the 
feedback and coupling sections on a subwavelength scale 

4>ij = ^oTij , (9) 



Chaos synchronization in networks of delay- coupled lasers 



6 



where Qq is the optical frequency of the laser. Since is larg^ one can for large delay 
consider the phases as parameters independent of the coupling delays. We thus choose 
all delays equal, but consider the phases as free parameters. 

One important peculiarity of coupled lasers is that the lasers may synchronize with 
a relative phase shift u 

E,{t)=e'^E,{t) . 

To explicitly treat this relative phase shift, we transform to the new variable E2 as 

E^it) = e'^E^it) . (10) 

After substituting this into Eqs. ([T]) and omitting the tildas for simplicity we arrive at 
the following rate equations for the fields 

E2 = ... + K22e'^'' E^it - t) + K2ie'^'^''^''^ Eiit - t) . (11) 

The artificially introduced parameter u helps in the analysis of synchronization, because 
we can discuss existence and stability of the synchronized solution Ei{t) = E2{t) in 
dependence on u, thereby treating synchronization of the lasers with a phase shift. 
Bringing Eq. (11) into the form of Eq. ([T| essentially yields the coupling matrix 



G 



(12) 



The row sum condition CI is then given by 

Kue"'"' + /€i2e^(*^^-") = 1^226'^''^ + ft2ie^('^^^+") . (13) 



Equation (13) can be interpreted as follows: If for a given set of coupling strengths and 



coupling phases there exists a phase shift u, such that Eq. (13) is fulfilled, then there 
exists an invariant synchronization manifold. The lasers can thus tune their relative 
phase shifts appropriately. The stability of the synchronized solution then determines 
whether synchronization will be observable or not. 

As a starting point for the stability analysis we consider the MSF for a laser network 
with row sum a. Since we focus on the large delay case (r = 1000 3> T^o) the rotational 
symmetry discussed above holds and the MSF A(re*'^) depends solely on r. Figure [2] 
depicts the MSF in the (cr, r)-plane for two different values of the pump current p = 0.1 
(panel (a)) and p = 1 (panel (b)). When the critical radius tq (solid line) lies below the 
diagonal line r = a (dotted line), the synchronized solution is chaotic for a network with 
this row sum, since A((t) > A(ro) = 0. When Tq = a, as occurs for instance in panel (b) 
for small values of a, the synchronized dynamics is periodic. In this case the solution 
A = corresponds to the Goldstone mode of the periodic orbit. 

I Typically, the optical frequency is of the order 10^'*s~^. We here use dimensionless units, where 
time is measured in units of the photon lifetime (e.g., Tp — 10~^^s). In these units typical optical 
frequencies correspond to a value of ~ lO'^. 
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Figure 2. MSF A(r) in the {a, r)-plane for p = 0.1 (panel (a)) and p — 1.0 (panel 
(b)). The solid line depicts the critical radius tq with A(ro) — 0. The dotted line is 
the diagonal line r — a. The blue dot corresponds to a parameter set that we will use 
in the following numerical investigations. 

a) Cj£r-"(D e) &^(2XP 

b) CjSr^©;) f) (D^Z^i) 

c) CjStzj© g) Cd) 

d) CjSrZjgQ 

Figure 3. Possible motifs for two delay-coupled lasers. Motifs a)-d) can (for certain 
parameters) exhibit zero-lag chaos synchronization and are discussed in more detail 
below. The motifs e)-g) cannot exhibit zero-lag chaos synchronization. 



In the following we will consider different coupling topologies of the two lasers. 
Figure |3] depicts all possible network motifs (up to exchange of 1 ^ 2) with more than 
one connection. The motifs on the left (a-d) can exhibit chaos synchronization in the 
limit of large delays, while those on the right (e-g) cannot |22] (trivially for motif (g)). 

We will now discuss the implications of the synchronization conditions C1-C3 for 
these motifs and compare the predictions with numerical simulations. We study two 
measures for synchronization. The first measure is the correlation coefficient of the laser 
intensities Ji, I2 

(ih - {h)) ih - {h))) 
P = — , (14) 

where (■) denotes the time average and (AJ^)^ denotes the variance of the respective 
intensity. Although, the correlation coefficient can in principle distinguish between 
identical synchronization (p = 1) and imperfect (generalized) correlation (p < 1), it is 
not the most sensitive measure for this purpose, because imperfect synchronization may 
still yield very large correlation p 1. 

To overcome this disadvantage of the correlation coefficient, we calculate the 
synchronization probability Ps- We define Ps as the probability that at any time t 
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the relative error between Ii{t) and l2{t) is smaller than a threshold e 

\hit)-h{t)\ 



We choose e 



0.01 in the following. 



< e 



(15) 



3.1. Motif (a) 

This case is the classical master slave configuration for chaos communication with lasers. 
It is also referred to as open-loop master slave configuration [35], since the receiver (2) 



G 



Kiie 



o«{<^21+m) 



has no self-feedback. The coupling matrix Eq. (|12|) is in this case given by 



K2ie-'-'^'^'''"-' 
The row sum condition CI then becomes 

This condition is satisfied if and only if Ku = k,21- The phase shift u can then compensate 
any choice of coupling phases (pu and 02i- The second eigenvalue of the coupling matrix 
is zero, such that the motif is optimal for chaos synchronization and condition C2 is 
always fulfilled. Whether the systems will synchronize or not then depends on whether 
A(0) is smaller or larger than zero, i.e., whether the MSF in Fig. [2] is negative (weak 
chaos) or positive (strong chaos) for r = 0. 

Note that since the coupling is unidirectional, the synchronization properties do 
not depend on the coupling delay. In particular, if the coupling delay is smaller than 
the self-feedback delay, anticipated synchronization [36] can be observed 



3.2. Motif (b) 

This setup consists of two unidirectionally coupled lasers with self-feedback and has been 
studied in different contexts. The importance of the coupling phases in this coupling 
scheme has been recognized in [20j. In this reference it was observed in an experiment 
that depending on the (relative) feedback phases the synchronization behavior ranges 
from perfect synchronization to an almost uncorrelated state. So far these experiments 
have not been sufficiently explained. It turns out that the experimental results can be 
well understood in the light of the synchronization conditions C1-C3: 
The coupling matrix corresponding to motif (b) is given by 





G 



and the row sum condition CI reads 



(16) 



Kue"^'' = K2ie"''^^^ ' + ^226 

Eliminating m, yields the following condition on the coupling strengths and phases 



«21 



\Kii - K22e 



"11 



+ K 



22 



2^11^22 COs($rel) 



(17) 
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Figure 4. (Motif (b) ) Synchronization probability Ps (left) and correlation coefficient 
p (right) as a function of (f>22 = ^tc\ and ki\. The dotted line corresponds to the 



synchronization condition CI (Eq. (17l). Other parameters: k\\ = 0.2, K22 = 0.09, 



0. 



with $ 



rel 



^22 



U 



. The phase shift u between the laser fields is then given by 



Arg (kii 



^^226 



where Arg denotes the complex argument. 

Condition C2 on the other hand reads for this case 

K22 < Ku = \cr\ , 

i.e., chaos synchronization is possible, if the self-coupling of laser 2 is weaker than that 
of laser 1. For a given set of parameters, synchronization is stable if condition C3 is 
fulfilled, i.e., if 



K22 < ro(/tii) . 

To illustrate these conditions, we consider the following parameter set 
0.2, K22 = 0.09, (i)ii=0, (i)2i=0 



Kn 



(19) 



(20) 



and vary ^^ei = 022 and k,2i. These coupling strengths (a, r) = (kh, K22) are marked 
by a blue dot in Fig. |2^ and are chosen such that synchronization is stable (condition 
C3 is fulfilled). Figure |4] depicts the correlation coefficient p and the synchronization 
probability Ps in the ($rei, K2i)-plane. The dotted line corresponds to Eq. (17), where 



condition CI is satisfied. This condition clearly coincides with high synchronization 
probabilities. 

As discussed above, high correlations are also possible without perfect 
synchronization. In the regions of high correlation and low synchronization probabilities, 
we observe generalized synchronization. Figure |5]^a) shows an example intensity time 
series observed in these regions. The parameters correspond to the blue dot in Fig |4j 
The blue dotted line (condition CI) in Fig |4] roughly marks the boundary between 
regions of high and low correlations. Below the line, we observe no synchronization 
at all. An examplary time series in this regime (corresponding to the white cross in 
Fig. |4|) is shown in Fig. |5|(b). The correlation in this case is very low. Interestingly, 
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200 400 600 800 1000 
(a) t 




200 400 600 800 1000 
(b) t 



Figure 5. Panel (a): Generalized synchronization for K21 = 0.15 022 = ^'roi = 
(blue dot in Fig. Panel (b): No synchronization for K21 — 0.15 (j)22 = $rci = f tt 
(white cross in Fig. |4]). Other parameters as in Fig. |4] 



the dynamics of the second laser has, in this strong high-frequency component. 

The occurrence of this high-frequency component can be understood as a result of the 
non-locking behavior. Since laser 2 does not lock to the signal of laser 1, the overall 
input signal of laser 2 is given by the interference of two (almost) independent chaotic 
signals, namely the signal from laser 1 and the self-feedback signal from laser 2. In this 
interference signal the high-frequency component is present due to a fast alternation of 
constructive and destructive interference. Indeed, by switching off the coupling {k,21 = 0) 
and calculating the intensity of an interference signal Jjntf = \Ei{t) + E2{t)\'^ one can 
already observe the high-frequency component. For non-zero coupling this interference 
signal drives laser 2, leading to even higher-order effects. 

In Fig. Qleft) there is a region on the dotted line ^rei ^ [0.27r, O.Stt], where the 
synchronization probability is low. In this region we observe mult ist ability between 
identical synchronization solution, a state of generalized synchronization similar to 
Fig. |5|^a), and an uncorrelated state similar to Fig. ^h). Which state is chosen depends 
in our deterministic simulations sensitively on the initial conditions. Including noise 
in the simulation results in spontaneous switching between the three states, albeit it is 
hard to distinguish between identical synchronization and generalized synchronization 
in the presence of noise. 

Thus although our synchronization condition C1-C3 give the correct existence and 
stability of the identical synchronized solution, we certainly cannot exclude the existence 
of other attractors. 

The experimental investigations in Ref. |i20j were performed under similar parameter 
conditions, i.e., K22 < (obeying condition C2). As $rei was varied, the correlation 
varied from almost perfect to almost no correlation. Varying $i.ei for a fixed value of 
K21 in Fig. |4] reproduces this behavior, provided the value of K21 intersects the dotted 
synchronization curve. Thus the experimental results can be understood as interference 
effects, which may or may not lead to the existence of a synchronization manifold, 
corresponding to high and low correlation, respectively. 
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3.3. Motif (c) 

For this motif the couphng matrix is given by 



G 



fi;2ie*(<^2i+n) Q 



i(0i2-«) 



and the row sum condition CI reads 



i{(j)l2-u) 



(21) 



(22) 



In an experiment with a simple bi-directional coupling, i.e., face-to-face coupling of the 
lasers, the coupling from laser 1 to laser 2 has the same coupling strengths as the reverse 
direction: 



Kl2 ~ /t21 • 



(23) 



We will first consider the more general case, which can be realized in an experiment 
using optical isolators and separate beam paths for the two coupling directions. We will 
then consider the natural condition Eq. (23) as a special case. 

Our aim now is to derive a condition that is equivalent to the existence of a phase 
shift u satisfying Eq. (22), i.e., we want to eliminate u from the equation. To simplify 
the discussion, we introduce two parameters: 

V = u-\ , and 9 = . (24) 

Note that f is a "free" parameter, since the phase shift u can be selected by the system. 
With these parameters Eq. (22) can be written as 



(25) 



For varying v the terms on the right hand side describe an ellipse in the complex plane 
with semi-minor axis a = \k,i2 — k,2i\ oriented along the real axis and semi-major axis 
b = |/ti2 + 't2i| oriented along the imaginary axis. In order for Eq. (25) to have a solution, 
the left hand side has to lie on this ellipse. Thus 



X 



ii<l>ii-e) 



y ■-- 



Re|^Kiie 
Iml'^iie^^'^" 



Kii COS(^</)ii - 

Kii sin(</)ii - 



have to obey 



62 



1 . 



(26) 
(27) 

(28) 



This equation is the desired synchronization condition corresponding to CI. It involves 



the coupling strengths Kh, Ki2, and K21 and the coupling phases 
can write it as an explicit equation for Kh 



>12' 



COS 



>ll 



sm 



^11 



L (K12 - K21, 



(Ki2 + K21) 



-1/2 



and (/)2i . We 



(29) 



Note that for the degenerate case /ti2 = K21, i-e., a = 0, the relevant condition becomes 



X 



0, 



and 



(30) 
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b = Ki2 + K21 
a = \ki2 — K2l\ 



X (Re) 



Figure 6. (Motif (c)) Ge- 
ometric visualization of the 
synchronization condition 



(Eq. (125|)). If the complex 
number Kne*^"^"^^^ lies on 
the ellipse, there exists an 
invariant synchronization 
manifold. 
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^ii/tt 
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Figure 7. (Motif (c)) Synchronization probability Ps (left) and correlation 
coefficient p (right) as a function of and kh. The dotted line corresponds to 



the synchronization condition CI (Eq. (29l). Other parameters: K21 = 0.2, K12 — 0.09, 

012 = 021 = 0. 



A necessary condition on the coupling strengths for the existence of a solution is 

— f^2l\ < f^U < fil2 + 1^21 ■ (31) 



To obtain the phase shift u, we split Eq. (25) into real and imaginary part and solve for 
V. This yields 

v = Argl + i ) . (32) 

\K21-I^12 1^21 + 1^12 J 

The phase shift u is obtained from the definition of v = u + (02i ~ 0i2)/2. 

We now discuss condition C2 and C3, which concern the eigenvalue of G other 
than the row sum. When G has a constant (complex) row sum a, this row sum is an 
eigenvalue of G. Since the determinant of a matrix is the product of its eigenvalues, we 
have the following equation for the second eigenvalue 71 of G 

|detG| = |/ti2/t2i| = |c"||7i| • (33) 
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On the other hand, there is only one entry ^216^^'^^'^^'^^ in the second row of G (see 



Eq. (21)), such that in the case of constant row sum we have 

|cr| = K21, and |7i| = K12 . (34) 
Thus chaos synchronization is possible if (C2) 

K12 < K21 = \cr\ . (35) 
For given parameters (given Tq) synchronization is stable if and only if (C3) 

K12 < ro(K2i) • (36) 
Again, we illustrate these conditions by numerical simulations. We choose 

Ki2 = 0.09, K21 = 0.2, 012 = 0, 021 = . (37) 

The coupling strength again corresponds to the blue dot in Fig. |2^, such that 
synchronization is stable if a synchronization manifold exists. Similarly to the case of 
motif (b), there are small regions of low synchronization probability on the dotted curve 
m Fig.|7|(left). These again correspond to regions of multistability and the simulations 
depend sensitively on initial conditions. 

Coming back to the most natural face-to- face coupling (Eq. ([23|)) with 

K12 = K21 (38) 



it becomes clear from Eq. (35) that chaos synchronization is not possible in this simple 



setup. To synchronize the two lasers in this motif one needs an asymmetric coupling. 
3.4. Motif (d) 

We now consider motif (d), where all possible couplings are present. The corresponding 
coupling matrix is given by 



G 



(39) 



K2ie'^'^^' ' K22e' 

The row sum condition CI is given by 

Kiie^"^" + /ti2e^(*i^-") = K2ie*('^2^+") + ^226^'^^^ . (40) 

Below we will discuss this condition in more detail. First, however, we turn to conditions 
C2 and C3 involving the eigenvalues a (row sum) and 7 of the matrix G. 



Assuming that the row sum condition Eq. (40) is satisfied, we solve Eq. (40) for 
^_i^^g«0ii a^nd replace this term in the coupling matrix (39). The resulting matrix then 
has the eigenvalues 

a = K22e'^^^ + K2ie^('^2i+'') and (41) 

7 = K22e''^^^ - /€i2e^('^^2-") . (42) 

The condition C2 (I7I < |cr|) gives (after some straightforward calculations) 

K12 — 2^12^22 COS(012 — 022 — u) 

< + 2K21K22 COS(021 - 022 + m) • (43) 
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The sufficient condition C3 similarly yields 

- 2K12K22 COS(012 - 022 - m) + I^l2 < ^o(o-) • (44) 

In contrast to the previous cases, it is not possible to eliminate the phase shift u between 
the lasers directly. We will, however, obtain a formula for u below. 



We now further discuss the row sum condition Eq. (40) and aim to eliminate the 
relative phase shift u. The calculation is similar to that done for motif (c), only more 
involved. Using the same parameters as for the case of motif (c) (Eq. ([24|)) 



V = , and 6 = , (45) 

we obtain 

ftue*^'^""^^ - K22e*('^22-^) = -Kise"™ + K2ie™ . (46) 

Again w is a free parameter, since it is proportional to the relative phase shift u, which 
can be selected by the system. This equation corresponds to the geometric problem 
visualized in Fig. |8| For varying v the terms on the right hand side describe an ellipse 
in the complex plane with semi-minor axis a = \ki2 — K21I oriented along the real axis 
and semi-major axis b = K12 + k,21 oriented along the imaginary axis. For the equation 
to have a solution the real and imaginary part of the left hand side have to lie on this 
ellipse. Thus 

X := Ku cos(0ii - 6) - K22 COS(022 - d) 

y := Ku sin(0ii - 9) ~ K22 sin(022 - 9) (47) 



have to obey 



-2+1^ = ^- (48) 



Equation (48) is the final condition, which has to be fulffiled in order for the 
synchronization manifold to be invariant. 

Movie 1 (animation.mpg) illustrates how the geometrical problem of Fig. |8] results 
in the dotted synchronization curves in Fig. |9} 

Figure [9] depicts the correlation coefficient and the synchronization probability in 
the (011, 022)-plane for a fixed set of coupling strengths and cross-coupling phases. On 



the dotted lines the synchronization condition CI (Eq. (48)) is fulffiled. The lasers 
exhibit strong synchronization on this curve. For this motif, the row sum a and the 
transversal eigenvalue 7 depend on the parameters 0ii and 022- Due to this dependence, 
the synchronization condition C3 is not satisfied everywhere on the dotted line and 
identical synchronization is only stable on parts of the curve. Additionally, we again 
have the effect of multistability as discussed for the last two motifs, which can result in 
low synchronization probability on the dotted curve. 

Similarly to the case of motif (c), the relative phase shift u can be found by solving 



Eq. (|46|) for v 

V = Are I 

'«21 - /tl2 '^21 + /tl 



Arg( +t ^ ) (49) 
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y (Im) A 



b = Ki2 + K21 
a = |ki2 — K21I 




Figure 8. (Motif (d)) Ge- 
ometric visualization of the 
synchronization condition 
(Eq. (|46|). If the sum of 
the two complex numbers 
Kne*^'^"-^) and K22e*("^^^"^) 
lies on the ellipse, there exists 
an invariant synchronization 
manifold. 
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Figure 9. (Motif (d)) Synchronization probability Ps (left) and correlation 
coefficient p (right) as a function of (pn and (/)22. The dotted line corresponds to 



the synchronization condition CI (Eq. (48 1). Other parameters: kh — 0.25, K12 — 0.1, 
K21 = 0.25, K22 — 0.15, (pi2 = 021 = 0. 



and using the definition oi v = u + (021 — (pu)/^. This phase shii^ can then be used in 
condition C2 and C3 (Eqs. (|43|) and (|34|)). 

In order for Eq. (48) to have a solution, the two vectors with respective lengths Kh 
and K,22 have to be able to reach the ellipse, i. e., the sum of the magnitudes Hu + k,22 
has to be larger or equal to the length of the semi-minor axis a = \ki2 — ^21 1- Similarly, 
the absolute value of the magnitude difference l^n — K22I has to be smaller or equal to 
the length of the semi-major axis b = K12 + k,21- We thus obtain two conditions for the 
existence of a solution 



1^22) 
1^22) 



> 
< 



^1^12 
[1^12 



1^21) 



(50) 



If and only if the coupling strengths fulfill Eqs. (50), there is a combination of phases 
such that condition (48) is satisfied. 



Chaos synchronization in networks of delay- coupled lasers 



16 



In many optical setups the forward and backward directions have approximately 
equal coupling strengths k,i2 ~ ^21- This holds for instance for a setup where the lasers 



are coupled via a common mirror (see Sec. 4.2) or where the lasers are coupled face-to- 



face. In this case the ellipse becomes a line along the y-axis stretching from —b to b. 
Thus Eq. (48) reduces to x = with the supplementary condition y'^ < b^. 



4. Laser networks 

We now move from the two-laser system to networks of ail-optically coupled lasers. 
The synchronization of laser networks is on one hand important for applications for 
instance in high power laser arrays, where the synchronization of optical phases yields 
an intensity I (x N'^ for interfering beams in contrast to the case of randomly distributed 
phases that gives I oc N ||351 [Ml- 

On the other hand, laser networks have been proposed as optical information 
processing systems HT] . Understanding stability properties of dynamical states in 
these networks is a necessary first step for controlling and utilizing these systems. 

For a network of ail-optically coupled lasers, the coupling terms in the Lang- 
Kobayashi rate equations are given by 

N 

E, = ... + Y,^kje'^'^Ej{t-T) . (51) 
i=i 

Here Kkj and (pkj are matrices describing the coupling strength and coupling phase, 
respectively, of the connection j — )■ k. Again the lasers may synchronize with relative 
phase shifts Uk, which we define with respect to laser 1 

Ei{t)=e'^'Ek{t) , (52) 



i.e., Ml = 0. Performing the corresponding transformation (as before in Eq. (11)) 

Ek{t) = e*"^Efc(t) (53) 
and omitting the tildas for simplicity we obtain the field equations 

N 

Ek = + e*('^'=^+"'=-"^)Ej- (t - r) . (54) 

j=i 

The necessary condition CI for synchronization is then stated as follows: There exists 
an invariant synchronization manifold if and only if there is a combination of relative 
phases Uk {k = 2, . . . , N), such that the complex row sum 

N 

ak = Y^ e^('^'=^+"'=-"^) (55) 

is independent of k. As we saw before, this condition is already very difficult to analyze 
for two lasers when all four connections are present. In an experiment it may be possible 
to actively control two or three phases accurately using for instance piezo positioning 
of mirrors, or passive wave-guides where the optical path-lengths can be controlled 



Chaos synchronization in networks of delay- coupled lasers 



17 



through an injection current |32]- However, it is virtually impossible to actively control 
more than a few phases in this way. For larger networks, optoelectronic |l3l HU H5] . 
incoherent optical feedback |^ Wf\ or coupling via a common relay [T3| [T6] may thus be 
a more promising coupling method. However, in certain cases all-optical coupling may 
be feasible in an experimental situation. We discuss such a feasible setup below. 

To illustrate the complexity of the synchronization condition CI in a simple case, 
we consider a system of bidirectionally coupled lasers. 



4.I. Rings of bi-directionally coupled lasers 

For rings of bi-directionally coupled lasers, the coupling matrix has the following 
principal structure 

Gi2 Gin 

G21 6*23 

G= G32 ■•■ ■•■ , (56) 

Gn-i,n 

Gni Gn,n-i 

with Gkj = t^kj 6**^'^'=^ "'"J-'. We then obtain the following row sums 

(72 



K12 e^(<^12+«l-«2) ^ f.^^ gi{<^iJV+«l-«iv) 



0"Af-l 



K^r , p*(</'iV-l,JV+"iV-l-W]v) I 1^ , i{(j)N-l,N-2+UN-l-UN-2) 

'^N-1,N ^ '^N-l,N-2(i ; 



ajv =kni e*(<^'^-i+"'^-"i) + kn n-i e^('^^>^-i+"'^-«^-i) 



(57) 



The necessary synchronization condition CI then corresponds to 2(A^ — 1) equations 
(each complex equation yields two real equations) 

cTi = 0-2, 0-2 = 0-3, ... aN-1 = ctn ■ (58) 

However, there are only — 1 relative phases uj {ui = 0) that the system can choose 
freely. Thus, in an experiment we need to control — 1 coupling phases to satisfy the 
synchronization condition CI. For a ring of = 3 elements, controlling two coupling 
phases may still be feasible, but this approach quickly fails for larger N. 



4-2. Coupling via a common mirror 

We will now discuss one promising all-optical setup that should, in principle, be robust 



to phase mismatches. Consider the setup [40j sketched in Fig. 10 The laser fields are 
coupled into a common fiber, which is terminated by a mirror. As before, we assume 
that the coupling delays are equal (the light paths are equally long on a cm scale), but 
allow for coupling phases, i.e., differences in optical path-lengths on wavelength scales. 
In this setup, the connection from each laser k to the mirror corresponds to a certain 
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mirror 



Figure 10. Coupling 
of multiple lasers via a 
common mirror. 



optical path length with a corresponding phase ijjk- The coupling phase from laser j to 
laser k is then given by 

(pkj = i'k + i^j ■ (59) 

Similarly, the connection from each laser to the mirror has an associated coupling 
strength Ck, which could in an experiment be controlled by an attenuator in the 
corresponding fiber. The coupling from laser j to laser k then has an effective coupling 
strength 

K-kj = Ck ■ Cj . (60) 
Under these conditions, the row sum is given by 

TV 

ak = Y^ e^('^^^+"^-"^) = Cfee^(^*+"'=) ^ e^^^^-^^^ . (61) 
i=i i=i 
The sum on the right hand side is a complex number independent of k. The row sum 
condition can thus be satisfied if the prefactor is also independent of fc, i.e., all coupling 
strengths are equal (cfc = c) and the relative phases are chosen by the system to giv^ 

^i=^k + Uk {k = 2,...,N) . (62) 

In essence, all lasers couple to the optical mean field and compensate for the difference 
in optical path length of their individual fiber by adjusting their relative phase shift. As 
long as the setup obeys Eqs. (59) and (60), the coupling phases tpk niay vary, e.g., due 
to thermal effects, and do not need to be controlled. 

Assuming that all coupling strengths are equal and the relative phases Uk are tuned 
appropriately by the system the coupling matrix is given by 

Gkj = e*('^''-+^^'+"'=-"^) = e^2^^- . (63) 

This matrix has one eigenvalue corresponding to the row sum 

N 

a = c2^e*2^^ (64) 



§ The special role of laser 1 stems from our choice ui = 0. 



Chaos synchronization in networks of delay- coupled lasers 



19 



and the — 1 transversal eigenvalues are zero 



In 







(n 



N -1) 



(65) 



such that this setup is optimal for synchronization. 

Recently, this coupling scheme has been proposed for optical information processing 
using multi-mode lasers |1D]. For multi-mode lasers, the same argument as above 
holds for each mode, such that the setup is, in principle, robust to phase mismatches. 
However, it may still be difficult to realize the assumptions Eq. (59) and Eq. (60) in an 
experimental setup. 



5. Conclusion and outlook 



We have discussed chaos synchronization conditions for ail-optically coupled lasers. In 
all-optical coupling, the coupling phases play a crucial role for the synchronizability. 
The condition of constant row sum corresponds to specific interference conditions, i.e., 
the input signals of each laser should interfere in such a way that each laser receives the 
same input signal, relative to its own phase. This corresponds to the existence of an 
identical synchronization manifold. Through interference, the phases may compensate 
for mismatches in the coupling strengths. 

Using these interference arguments we have explained experimental findings |20j and 
discussed necessary and sufficient conditions for synchronization of all network motifs 
which contain two ail-optically coupled lasers. 

Further, we have considered synchronization of larger laser networks, and singled 
out the difficulties that arise in all-optical coupling schemes due to the interference 
conditions. We predict that a setup of all-to-all coupling via a common mirror may under 
certain conditions be robust to phase mismatches and thus be optimal with respect to 
stability of the synchronized chaotic dynamics. 

In a broader context, these results might also be relevant for other networks where 
phase-sensitive couplings play a role, e.g., networks of Stuart-Landau oscillators with 
complex coupling constants a = Ke^^ j7]. These are generic models representative for 
a large class of oscillator networks. 
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